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Abstract 

We investigate the 5D massless fermionic fields within the standing wave braneworld 
model. We show that in the case of increasing warp factor there exist localized left spinor 
field zero modes on the brane, while right fermion wave functions are not normalizable. 



PACS numbers: 04.50.-h, ll.25.-w, 11.27.+d 



The brane models [HE] has attracted a lot of interest recently with the aim of solving 
several open questions in modern physics. A key requirement for realizing the braneworld 
idea is that the various bulk fields be localized on the brane. For reasons of economy and 
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avoidance of charge universality obstruction [3] one would like to have a universal gravitational 
trapping mechanism for all fields. However, there are difficulties to realize such mechanism 
with exponentially warped space-times. In the existing (l+4)-dimensional models spin and 
spin 2 fields can be localized on the brane with the decreasing warp factor [2], spin 1/2 field 
can be localized with the increasing factor [4], and spin 1 fields are not localized at all [5]. For 
the case of (l+5)-dimensions it was found that spin 0, spin 1 and spin 2 fields are localized 
on the brane with the decreasing warp factor and spin 1/2 fields again are localized with the 
increasing factor [B] . There exist also 6D models with non-exponential warp factors that provide 
gravitational localization of all kind of bulk fields on the brane [TJ, however, these models require 
introduction of unnatural sources. 

To solve the localization problem recently we had proposed the standing wave braneworld 
model [S] , which is generated by collective oscillations of gravitational and scalar phantom-like 
fields (similar to [5]) in 5D bulk space-time. The metric of the model in the case of increasing 
warp factor has the form: 

d s 2 = e 2a|r| (dt 2 - e u dx 2 - e u dy 2 - e~ 2u dz 2 ) - dr 2 . (1) 

Here a = a/A/6 > 0, where A is 5D cosmological constant, and 

u{t,r) = B sm{ut) £(r) , £(r) = e" 2a|r| y 2 (^e" a|r| ) , (2) 

where B is a constant, u denotes the oscillation frequency of the standing wave, and Y 2 is 
second order Bessel function of the second kind. The solution ([1]) describes the brane at r = 0, 
which undergoes anisotropic oscillations and sends waves into the 'sea' of phantom-like scalar 
field in the bulk. 

The r-dependent factor £(r) of the metric function ([2]) has finite number of zeros along 
the extra coordinate axis, which form the nodes of the standing wave. These nodes can be 
considered as 4D space-time 'islands', where the matter particles can be bound [8]. One of 
nodes of the standing wave must be located at the position of the brane. This can be achieved 
by imposing the condition 

Y 2 Q = (3) 

fixing the value of the frequency of standing wave, u, in terms of the curvature scale a. In what 
follows we assume: 

- = Z x ~ 3.3842, (4) 
a 

where Z n denotes the n-th zero of the function Y 2 . 

The standing wave can provide localization of the matter particles with energies much 
smaller than u. One of the classical analogs of this mechanism is the so called optical lattice 
with standing electromagnetic waves [10] . In the previous papers [TT] we had demonstrated the 
explicit localization of scalar, vector and tensor fields zero modes in standing wave braneworld 
with increasing warp factor. In this article we investigate the localization problem for massless 
fermions within this model. 
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For Minkowskian 4x4 gamma matrices ({7°, 7^} = 27/ Q/3 , Greek indices, a, /?,... = £, cc, y, z, 
numerate 4D coordinates) we use the Weyl basis, 

1' = ( ? ) ■ -'•=("' ^' ) • ^ = *»W = ( -°, ) • < 5 > 

where I and a* (z = x, y, z) denote the standard 2x2 unit and Pauli matrices respectively. 

Let us recall that four-component columns represent fermions in 5D, and that 5D gamma 
matrices can be chosen as: 

Y A = h A T A , = (7 4 ,7 x ,7 y ,7 2 ^7 5 ) , {T A ,T B } = 2g AB , (6) 

where capital Latin indices, A, B,... = t,x,y, z,r, stand for 5D space-time coordinates and 
A,B, refer to 5D local Lorentz (tangent) frame. So according to (JHD the curved-space 
gamma matrices T A are related to Minkowskian ones (jSJ) as: 



(7) 



The fiinfbein for our metric ([T]) , 

h A = (e a ' r ' e a ^ +u ^ 2 e a ' r ' +U//2 e a ^~ u l) 

h AA = g AB h A = ^ e -a\r\^_ e -a\r\-u/2^_ e -a\r\-u/2^_ e -a\r\+u^_ 1 ^ (g) 

h\ = r] AB h BA = (e- alrl ,e- alrl - u/2 ,e- alrl - u/2 ,e~ alrl+u ,l) , 
hiA = VaS^a = (e a l r l,-e a l r l +M / 2 ,-e a l r l +u / 2 ,-e a l r l- u ,-l) , 

is introduced through the conventional definition: 

9ab = VABhfhl . (9) 

The 5D Dirac action for free massless fermions can be written as: 

S = J d 5 Xy/g~ W (x A ) r M D M V (x A ) , (10) 
where the determinant yfg of the metric ([1]) is: 

Vg = e 4alrl i (11) 

and covariant derivatives are defined as follows: 

D A = d A + -nfTzT c . (12) 



In the last expression fij^ N denote the spin-connections: 



h NM (d M h N N - d N hl\ - h NN d M h% - d N hfA - 



r\MN _ n NM _ 1 

\L M — —\L M — — yi yUM^N ~ U NII>M J ~ n y u M IL N ~ u N> l M 

-hih™h?» (d P h QA - d Q h PA )] . (13) 

The non-vanishing components of the spin-connection in the background ([I]) are: 

nf = - (e^V , 

of = nf = - ( e alrl+u/2 )' , 

Of = - (e a|r| - n )' , (14) 
- d (e u/2 ) 

r\xt _ r\yt _ \ ) 

x ~ y ~ dt ' 



n 



zt 



d (e~ u ) 



where primes denote derivatives with respect to the extra coordinate r. 
The corresponding to ( TTUj) 5D Dirac equation reads: 

iT A D A q = i (VD^ + T r D r ) 9 = 0. (15) 

For the bulk fermion field wave function we use the chiral decomposition: 

{x\ r) = ipL (x v ) A(r) + tfj R (x u ) p(r) , (16) 

where A(r) and p(r) are extra dimension factors of the left and right fermion wave functions 
respectively. We assume that 4D left and right Dirac spinors, 

correspond to zero mode wave functions, i.e. they satisfy free Dirac equations: 

i-fd^L = i-fdpipn = . (18) 

Apart from the massless states ipi and i/jr, the 5D Dirac equation also have solutions cor- 
responding to massive fermions. In the single brane models the masses of the bounded massive 
states are typically of order of the energy scale a, characterizing the brane as a topological 
defect in higher-dimensional space-time. These states are very heavy and we do not consider 
them here. 

The solutions of (fl8j) in our representation ([5]) can be written in the form: 

lj) R ( X U ) = ( f J e -i{Et- Px x-p y y- Pz z) ? 





c 

L 



IpL (X U ) = ( r ) e -^t- Px x-p y y-p z z) ^ ^ 
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where the constant 2-spinors L and R satisfy the relations: 

(E + (7%) L=(E- cr l Pl ) R = . (20) 

When the frequency u of standing waves in the background metric (OQ) is much larger than 
frequencies associated with the energies E of the fermions on the brane, 

u > E , (21) 

we can perform time averaging of the oscillatory functions in the Dirac equation ( Tl5|) . Explicit 
expressions for time averages where found in [TT] : 

<«> = w) = (§) = (^p) = > ( ebu ) = J o mm) , (22) 

where b is a constant and Jq is the modified Bessel function of zero order. Time averages of the 
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(23) 



where sgn(r) is the sign function. 

Now the equation (fl"5]) can be written in the form: 

% \fdt + <e M/2 ) (7*9, + fdy) + (e~ u ) 7 Z <9 Z ] * = e a|r| 7 5 [2a ^n(r) + d r ] * . (24) 

Using the solutions of free equations (fl9l) and the relations (|20|) it can be rewritten as the 
system: 

f -e^[2asgn(r) + d r ] tr'P^r) \/ p(r)fl \ , , 

V -^Pi(r) e a l r l [2a s#n(r) + <9 r ] ^ V A ( r ) L / ' 

where we have introduced functions P%{r): 

P x {r) = {(e-^ 2 )-l) Px =[I (\B\t(r)/2)-l)p x , 
Py{r) = ((e^ 2 )-l)p y =[I (m(r)/2)-l]p y , 
Pz(r) = ((e")-l)p z = [I Q (\B\t(r))-i\p t , (26) 



P 2 (r) = P 2 + P 2 + P 2 
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These functions can be considered as the components of P(r): 

P 2 {r) = P 2 x +P 2 + P 2 , 

which we call V- dependent momentum'. 

From the second equation of the system f[2"5]) it is straightforward to find 

p{r)R = e^^^ [ 2a s 9<r) + d r ] \{r)L . 



(27) 



(28) 



Inserting ([281 into the first equation of (|25|) and multiplying the result by a 1 Pi, we receive the 
second order differential equation for the function A(r): 



5a sgn(r) — — 



A' + 



P' 



4a5(r) + 6a 2 - 2a sgn{r) p 2 e " 2a|r| 



P 



A = . 



(29) 



Now let us investigate this equation in the two limiting regions: far from and close to the 
brane. 

Close to the brane, r — > ±0, the V -dependent momentum' ( )26|) behaves as: 



P(r) 



Ar 2 + 0{r 3 ) 



where A is constant, and the equation ( 1291) takes the following asymptotic form: 



A" + 



5a sgn(r) 



4a 

4a5(r) + 6a 2 sgn{r) 



A = . 



This equation has the unique nontrivial solution: 

A(r)U ±0 = Ce- 2a H 



(30) 



(31) 



(32) 



where C is a constant. 

As it follows from (128]) and ( 132]) . in our setup the right fermionic modes are absent on the 
brane: 

p(r)U ±0 = 0. (33) 

Such different behavior of the left and right massless fermions on the brane is not surprising, 
since in our model the effective mass term in (1241 is of 7 5 -type, 

(34) 



m{r) = 2a^f 5 sgn(r)e a ^ 



with the gap: 

|m(r) - m(— r)| = 4a7 5 e a|r| . (35) 

In the second limiting region r — > ±oo the function P'/P vanishes and the equation ( 129]) 
obtains the asymptotic form: 



A" + 5a s#n(r)A' + 6a 2 A = 



(36) 
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with the solution: 

A(r)U ±0O ~e- 3a M . (37) 

Using ( 1371) from the relation ( 12 8j) we find the asymptotic behavior of the extra factor of the 
right fermion wave function: 

p(r)U ±00 ~e- 2a l r l . (38) 

The condition for the localization of a field on the brane is that its wave function in the 
extra dimensions be normalizable, or that its action integral over r be convergent. We had 
found that in our model the extra dimension part A(r) of the bulk left spinor wave function 
(fT6j) has the maximum at the origin, 

A(r)| r=0 = C, (39) 

falls off from the brane, and turns into the asymptotic form (1371) at the infinity. When r — > 
oo the determinant (fTTj) in the action integral for 5D spinor fields (fTUl) increases as e 4a ' r L 
However, extra dimension factor of left fermions according to (l3Tj) contribute e _6a ' r ' and overall 
r-depended part in ffTUl) decreases as e _2a ' r L So in the case of left fermions the integral over r 
in (fTUj) is convergent, i.e. zero modes of left fermions are localized on the brane. 

At the same time according to (|33|) right spinor zero modes does not exist on the brane and 
due to (138]) at the infinity the extra dimension part of right fermions p(r) decreases as e~ 2a ' r '. It 
means that for right fermions integral over r in (TTOT) diverges and zero modes of right fermions 
are not normalizable on the brane. 

To conclude in this latter we have studied localization problem of the massless fermions in 
the 5D standing wave braneworld model. We have found that in the case of increasing warp 
factor there exist left spinor field zero modes localized on the brane, while right fermionic wave 
functions are not normalizable. 
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